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ABSTRACT 

We use data from the first 100 deg^ field observed by the South Pole Telescope (SPT) in 2008 to 
measure the angular power spectrum of temperature anisotropics contributed by the background of 
dusty star-forming galaxies (DSFGs) at millimeter wavelengths. From the auto and cross-correlation 
of 150 and 220 GHz SPT maps, we significantly detect both Poisson distributed and, for the first 
time at millimeter wavelengths, clustered components of power from a background of DSFGs. The 
spectral indices of the Poisson and clustered components are found to be Q:f5o-220 ^ ^-^^ ± 0.23 and 
0^150-220 — 3.8 ± 1.3, implying a steep scaling of the dust emissivity index (3^2. The Poisson and 
clustered power detected in SPT, BLAST (at 600, 860, and 1200 GHz), and Spitzer (1900 GHz) data 
can be understood in the context of a simple model in which all galaxies have the same graybody 
spectrum with dust emissivity index of /3 = 2 and dust temperature Td = 34 K. In this model, half of 
the 150 GHz background light comes from rcdshifts greater than 3.2. We also use the SPT data to 
place an upper limit on the amplitude of the kinetic Sunyaev-Zel'dovich power spectrum at £ = 3000 
of 13 /xK^ at 95% confidence. 

Subject headings: cosmic background radiation — galaxies: abundances — large-scale structure of 
universe — submillimcter: diffuse background — submillimeter: galaxies 
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1. INTRODUCTION 

The cosmic infrared background (GIB) is produced by 
thermal emission fr om dust in galaxies o v er a very broad 
range in redshift (iLagache et al.l l2005t iMarsden et al.l 
l2009t) . The dust grains, ranging in size from a few 
molecules to 0.1 mm, absorb light at wavelengths smaller 
than their size, and re-radiate it at longer wavelengths. 
Sufficient absorption occurs to account for roughly equal 
amounts of energy in the GIB and in the unprocessed 
starlight that inakes up the optical/UV background 
()Dwek fc Arendt|[l998l : iFixsen et al.1[l998l ). 

The spectral shape of the graybody emission from the 
dust, steeply rising with increasing frequency to a peak 
near ^ 3000 GHz (100 /im), results in a very peculiar 
luminosity- redshif t relation for the g alaxies contributing 
to the GIB (e.g., M. 

ain et all I2002D . For example, an 
object with emission that peaks at a typical rest-frame 
frequency of 3000 GHz, has an observed v = 220 GHz 
(A = 1.4 mm) flux that is nearly independent of red- 
shift for 1 < z < 13, as the increased luminosity due to 
increased rest-frame frequency compensates for the cos- 
mological dimming with redshift. Since we do not expect 
significant emission from redshifts anywhere near z > 13, 
measurements of the GIB at 220 GHz arc sensitive to 
the complete history of emission from dusty star-forming 
galaxies (DSFGs). 

Berkeley Center for Cosmological Physics, Department of 
Physics, University of California, and Lawrence Berkeley National 
Laboratory, Berkeley, CA 94720, USA 
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iPuget et all ()1996D detected the CIB for the first 
time with data from the Far Infrared Absolute Spec- 
trophotometer (FIRAS) on the Cosmic Background Ex- 
plorer (COBE) satellite. The CIB was first par- 
tially resolved into light from individual galaxies in 
350 GHz (850 fim) surveys using the S CUBA camera on 
the James Clerk Maxwell Telescope (iSmail et al.l 119971: 
iHughes et all [19981 lEales et alj [19991 ) . iChapiiiaii et alj 
( 200l " found that sources with 350 GHz fluxes in the 2 - 
15mJy range have a mean redshift of z ~ 2.2 and a 
maximum of z = 3.6. 

Progress is being made in understanding the nature 
of the sources o f the CIB. Using deep Spitzer data, 
iDole et all (f2006l) stacked 4.3 THz (70 /^m) and 1.9 THz 
(160 /im) flux from more than 19000 12.5 THz (24 
/xm) sources to determine that these sources contribute 
more than 70% of t he CIB at 4.3 THz and 1.9 THz. 
IMarsden et al.l (j2009() perform a similar analysis with 
Balloon-borne Large-Aperture Submillimeter Telescope 
(BLAST) maps to show that Far-Infrared Deep Ex- 
tragalactic Legacy (FIDEL) catalog sources are con- 
tributing ~ 95% of the FI RAS-determined backgr ound 
flux at 600 GHz (500 /im). iLe Floc'h etHI (|2005l ) and 
ILagache et al.l (j2005f ) show that the fraction of the CIB 
originating from sour ces at high red s hift d ecreases with 
increasing frequency. IMarsden et al.l (I2009D confirm this 
and find that 60% of the CIB originates from z > 1.2 
sources at 600 GHz. By including sources detected with 
the assistance of gravitational lensing magnification, it is 
possible to account for close to 75% of the 350G Hz back- 
ground (jEMneFaD [Till; [KWdiiniFin [lOOl^ Signif- 
icant advances in identifying the sources responsible for 
the CIB are expected with Herschel/SVlRWs 600, 860, 
and 1200 GHz channels and their resulting catalogs. The 
full-sky surveys of the Planck/Y{F1 will complement Her- 
schel by fully exploring the bright and rare end of the 
flux distribution. 

Observations of the diffuse background complement 
those of individual sources for two reasons. First, due 
to limitations on the resolution of current instruments, 
the bulk of the light at frequencies lower than 12.5 THz 
is emitted by sources at or below the confusion limit. 
Second, the clustering properties of the CIB are sensitive 
to the relation between dust-enshrouded star-forming 
galaxies and the dark-m atter halos in which they reside 
(jHaiman fc Kno^j 120001 ) . This connection is critical to 
any theoretical understanding of the history of star for- 
mation. Clustering information can also come from very 
large source catalogs, but the number of sources required 
to compete with di ffuse clustering me asurements is in the 
tens of thousands (jKnox et al.ll200l"[ ). 

Predictions for the clustered CI B fluctuatio n s date 
bac k to the p i oneeri ng work of IBond et "all (|1986f ) 
and iBond eFaLl (|1991D . Interest was revived with the 
SCUBA measurements, the discovery and characteriza- 
tion of the background from COBE data (jPuget et al.l 
Il996t iFixsen et allll998[ ). and the prospect of far-IR to 
millimeter wavelength sky maps to come from a number 
of cosmic microwave back ground (CMB) mapping 
proje cts, including Planck ijThe Planck Collaborationl 
|2006( ). Predictions for the clustered and Poisson power of 

Usually defined as the flux at which sourc es reach a surface 
density of 1 object per 40 beams IICondonlll974IV 



"unresolved SCU BA sources" have been made by a num 
ber of authors ([Scott fc Whitd Il999t iHaiman fc Kno^ 
2000t iKnox et al.l 120011: iMagliocchetti et al.l 12001 



Perrotta et all 120031: ISong et al.l 120031: INegrello et al l 



20071: iFerna^^ -Conde et al. 1 120081 : IRighietal.ll2008D . 

Models generically predict a clustering amplitude on 
the order of a few percent of the intensity of the mean 
background and angular scales of a few arcminutes, 
which can be understood by the following argument. A 
few arcminutes corresponds to several comoving Mpc at 
cosmological distances. The variance in density fluctua- 
tions is of order unity on Mpc scales, and the fluctuations 
in the galaxy distribution will be larger by some bias fac- 
tor. The integrated fluctuations along the line of sight 
probe of order a few hundred independent cells; the fluc- 
tuations are thus suppressed by the square root of a few 
hundred. This leads to an expected clustering amplitude 
on the order of a few percent of the mean. 

The first detection of diffuse infrared background 
clustered power came from 160 nm. Sp itzer data 
(jGrossan fc SmootI 120071 : ILagache et al.ll2007[) . This was 
followed recently by BLAST which measured the clus- 
tered power near the energetic peak of the CIB at 600, 
860, and 1200 GHz (500, 350, and 250 ^im) (jViero et all 
l2009f ). Here, we report on the first detection of CIB 
clustered power at millimeter wavelengths, and compare 
models of the CIB with measurements of this clustered 
power over more than a decade in frequency. 

This paper is one of a set of companion papers present- 
ing early results from South Pole Telesco pe (SPT) obser- 
vati qns that include iVieira et al.l (j2009l hereafter V09) 
and iLueker et all (|2009l . hereafter LIO). V09 present 
detections and spectral index measurements of point 
sources that are sufficiently bright to be resolved from 
the diffuse background. LIO present fine scale CMB 
anisotropy bandpowcrs measured with SPT data and re- 
port the detection of the thermal Sunyaev-Zel'dovich 
(tSZ) power spectrum in a "DSFG-subtracted" linear 
combination of the 150 and 220 GHz maps. In support of 
the interpretation of the LIO power spectrum, this paper 
constrains the amount of residual CIB clustered power 
that can remain in that map despite the DSFG subtrac- 
tion. 

To interpret our clustering signal as arising from the 
clustering of DSFGs, we need to consider other possible 
sources of fluctuation power as well. We address the 
impact on our results of the predicted Poisson power from 
radio galaxies. We also use a number of lines of argument 
to conclude that the contributions from the clustering of 
radio galaxies are negligible for both the 150 and 220 GHz 
maps. 

A difficult contaminant for the interpretation of our 
data is the kinetic Sunyaev-Zel'dovich (kSZ) effect, which 
has a highly uncertain amplitude. We discuss implica- 
tions of the kSZ component for our determination of the 
clustering signal, and also provide an upper limit on kSZ 
power sufficiently tight to rule out the highest predictions 
in the literature. 

The paper is organized as follows. We describe the ex- 
ternal galaxy emission models we use in Section 2. In Sec- 
tion 3, we describe our modeling of other signals, namely, 
the CMB, tSZ, kSZ, and Galactic cirrus emission. In Sec- 
tion 4, we describe the data and its processing to 150 GHz 
auto spectrum bandpowers (referred to by 150 x 150), 
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the 150 and 220 GHz cross-spectrum (150 x 220), and 
220 GHz auto spectrum bandpowers (220 x 220). In Sec- 
tion 5, we fit the data to the power spectrum templates 
described in Sections 2 and 3. In Section 6, we discuss 
the imphcations of these results when combined with 
those of other experiments and compare the results with 
models for the properties of the sources composing the 
millimeter-wavelength to far-infrared background light. 
In Section 7, we calculate the expected level of resid- 
ual contamination from DSFGs in the DSFG-subtracted 
map of LIO. We summarize the main results in Section 
8. 

2. MODELING THE ANGULAR POWER SPECTRA OF 
LIGHT FROM EXTERNAL GALAXIES 

The bandpowers in LIO are reported as I?^, where 



In this work, we present constraints on measured powers 
in terms of either and T>i depending on the source. Cg 
is constant for a Poisson distribution of sources, and is 
well suited to describing the Poisson point source power. 
On the other hand, T>i is nearly constant for the kSZ, 
tSZ, and clustered point source components at the angu- 
lar scales of interest in this work. We generically refer to 
the tSZ, kSZ, and clustered contributions to the power 
spectra as "flat" components due to their weak depen- 
dence of "Di on t — weak compared to both the CMB tem- 
perature fluctuation power spectrum in the damping tail 
region and the Poisson power spectra. 

2.1. Poisson 

The discrete nature of galaxies leads to fluctuations 
in their number in any comoving volume that are often 
approximated as following a Poisson distribution. As- 
suming this distribution, the resulting power spectrum 
depends on the flux distribution function as follows: 

/■■Scut J AT 

Ci = S^—d{lnS), (2) 

and is independent of £. Here {dN/dS)AS is the number 
of sources per unit angular area in a flux bin of width 
AS-. 

In Figure [l] we plot the integrand in Equation [2] for 
models and for the measured source counts in V09. Note 
that with the logarithmic a;-axis, the Poisson power is 
proportional to the area under the curve. By "spec- 
tral index" in this work we mean a^-^-^^ the rela- 
tion S'l/j = 5't/3 (1/1/^^2)"""^""^ ■ The DSFG source -count 
model predictions are from iNegrello et al.l ()2007[) ^° ex- 
trapolated to SPT frequencies with two broken power 
laws since the low-z population is expected to have a 
spectral index that steepens with increasing frequency 
and the high-z population is expected to have a spec- 
tral index that decreases with increasing frequency (see 
Figure [7]). For the high-z (low-z) populations, we use 
"220-350 = 3.2(3.0) and ai5o-220 = 3.5 (2.0), with the 

These source counts are based on the physical model of 
IGranato et al.l I I2004I ) for high-z SCUBA-like sources (that we call 
high-z DSFGs in Figure [TJ and on a more phenomenological ap- 
proach for late-type starburst plus normal spiral galaxies (that we 
call low-z ULIRGs in Figure [TJ. 



shallower indices at low redshift due to synchrotron and 
free-free emiss i on. The radio galaxy source counts of 
Ide Zotti et al.l ()2005[ ) are also shown, extrapolated from 
150 to 220 GHz with ai5o-220 = -0.5. The vertical 
dashed line in the left panel indicates the flux thresh- 
old used at 150 GHz to determine the point source mask 
that is used in both the 150 and 220 GHz SPT maps. 
Also shown are source counts from SPT data (V09) sep- 
arated into DSFG and AGN populations. 

According to the models, radio galaxies and DSFGs 
have very different source-count flux distributions. While 
the radio galaxies dominate the SPT detections, their 
contribution to the Poisson power (in maps with detected 
sources masked) is subdominant to those from the DS- 
FGs. Below S ^ 3 mJy at 150 GHz, we see contributions 
from a large population of DSFGs extending in redshift 
from z 1 to much higher redshifts. We expand our con- 
sideration of radio galaxy contributions in Section [2.31 

2.2. Clustering of dusty star-forming galaxies 

In this section, we calculate the clustering properties 
of DSFGs from models and compute the resulting power 
spectra. We consi der two models in thi s work, the rich 
empirical model of iLagache et al.l (|2003| ). as updated in 
iLagache et al.l ()2004[) ^^. which we henceforth refer to as 
the LDP model, and a more simplistic "single spectral 
energy distribution (SED)" model. The advantage of 
the single-SED model is that it has a limited set of pa- 
rameters which can be varied to match both the SPT 
and external data. The template shapes that are used 
in the analysis of the SPT and BLAST band power 
data are presented in Section [5l Examples of other ef- 
forts to model the power from clustered poin t sources in- 
cludelSeheal et al.l (|2010l hereafter S09) and lRighi et"all 
(I2OO8II 

Foll owing iHaiman fc Knoxl (|2000[ ) and iKnox et al] 
()2001|) , we calculate the clustering properties of the back- 
ground light by assuming that the emissivity density, j, 
is a_biased tracer of the dark-matter density fluctuations, 
5j/j = bSp/p. In a flat universe, the intensity, I^, is re- 
lated to the comoving emissivity density via 

/.= /rfzf^a,>,z)(;i + ^), (3) 
J dz V Ji^^z) J 

where x is the comoving distance to redshift z and a — 
1/ (1 -I- z) is the scale factor. 
For an angular power spectrum deflncd via 

{5Iv,lm5Iv' ,l'm') = Clj^,iSii>Smm', (4) 

we get 

CL',i^ j ^'^^a'b^3{y.z)-j{v',z)PM{krz)\k^iU. (5) 

Here, we have assumed a spatially constant and 
frequency-independent bias, b. We further assume that 

PM(fc) = Po{k)G\z), (6) 

In order to take into account constraints on number counts at 
1.9, 4.3, and 12.5 THz (160, 70, and 2i^lm) obtaine d from Spitzer 
observations llPapovich et al.ll200l : IDole et al.ll2004l) 




Fig. 1.— S^dN/dS vs. flux (S) at 150 GHz (left panel) and 220 GHz {right panel). Source counts of iNeerello et all l|2007D 
from high-2 {z > 1) DSFGs {blue dotted curve), low-2 DSFGs {purple dashed curve, often caUed ultra-lumi nous infrared galaxies 
(ULIRGs) due to their detection at IR wavelengths) and the sum of these two {black solid curve). The Ide Zotti et al.l (|2005D 
model radio galaxy curve is also plotted {red dot-dashed curve), extrapolated from 150 GHz with a = —0.5. We extrapolate 
the Negrello model predictions as described in the text. The vertical dashed line indicates the flux threshold at 150 GHz set for 
masking of point sources from both maps. At the bright end, the source counts of V09 due to radio galaxies {filled triangles) 
and DSFGs {filled circles) are also plotted. 



where G{z) is the linear-theory density contrast growth 
factor normalized to unity today and Po{k) is the power 
spectrum today. 

Nonlinear effects (which wc do not include) are ex- 
pected to alter the shape and amplitude of the power 
spectrum on the relevant scales by boosting power on 
smaller angular scales. We have ignored these corrections 
due to our limited sensitivity to the shape of the clus- 
tered power spectrum and its model dependence. How- 
ever, the amplitude correction is important to keep in 
mind when interpreti ng our results. In the halo model of 
large-scale structure (jSeliak et al.ll200l[ ) , the power spec- 
trum of galaxies is the sum of a one-halo term and a two- 
halo term, depending on whether it arises from a pair of 
galaxies residing in one halo or two different halos. The 
extra power at small angular scales, above expectations 
from linear perturbation theory, arises from the one-halo 
term. The bias factors of the halos in which the galax- 
ies reside scale the amplitude of the two-halo term, and 
usually an effective bias is defined as a weighted aver- 
age over the bias factors of all the contributing halos. 
The bias wc define above is, therefore, different than a 
halo-model bias. Fo r example, in the a nalysis of 860 GHz 
maps f rom BLAST.IViero et al.l (|2009h find 6 = 3.9 ± 0.6 
for the lLagache et al.l ()2004D niiodel using a linear-theory 
calculation, and an effective bias of 6 = 2.4 ± 0.2, us- 
ing a halo model. Curiously, the shape of their power 
spectrum agrees better with the linear-theory shape than 
with any halo-model shape they manage to derive from 
their model. 

We express sky brightness fluctuations in terms of the 
departure from the mean CMB temperature that would 
give the equivalent change in brightness. The power spec- 
tra are related via 

/ dB \ / dB 



where B^, is the Planck function, the derivatives of which 
are evaluated at T = Tcmb- 

The emissivity density is related to differential source 
counts via 

= + 'iSS^^ (8) 

= (l + z) 

/.M^)5(x,.)e(5-..„o^^(|)(9) 

where Q{x — y) is the Heaviside function equal to 1 for 
X < y and for a; > y. We take (PN/{dln{L) dz) 
from the empirical LDP source-count model, designed 
to accommodate all observational constraints on source 
counts and Poisson fluctuation power from the mid-IR 
to submillimeter wavelengths as economically as possi- 
ble. In this model, there are two different spectral types: 
"normal" and "starburst" galaxies. The latter undergo a 
rapid evolution of the luminosity function between z = 
and z ^ I and provide the large majority of the back- 
ground light. 

In order to explore the broader implications of the 
SPT constraints on DSFG clustered power, wc use an 
economical empirical model of the mean comoving emis- 
sivity density as a function of redshift. We use this 
"single-SED model" rather than the LDP model since 
the LDP source counts are only publicly available for a 
fixed parameter set. In our single-SED model, we assume 
that all galaxies have the same grcybody SED described 
by fu oc By [Tj) across the far-infrared to millime- 
ter wavelengths. Here, /3 is the emissivity spectral index 
of thermal graybody dust, B^ is the Planck function, 
and Td is the effective dust temperature. We replace the 
mid-infrared exponential decline on the Wien side of the 
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graybody with a power-law decline oc jy^^mid-m i^y 
matching the two functions with a smooth gradient at the 
frequency v' that satisf i es d In f,y /dhiv' = — amid-iR, as 
is done in IB lain et al.l (|2003f ). This replacement phe- 
nomenologically accounts for emission from dust with 
temperature greater than T^, which becomes important 
beyond the peak of the emission from the dust with tem- 
perature Td. In Figure [2j we plot this SED for a specific 
choice of parameter values and show the range of rest 
frequencies probed by observations at SPT, BLAST, and 
Spitzer frequencies. 
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Fig. 2. — Assumed spectral shape of galaxies for our singie- 
SED model of the millimeter- wave background {black dotted), 
with our fiducial parameter values Td = 34 K, /3 = 2, and 
Qmid-iR = 2. Curves indicate the range of rest-frame fre- 
quencies of light emitted at redshifts between 1 and 5, for the 
observed frequencies of 1900 GHz {triple-dot-dashed green), 
600 GHz {solid black), 220 GHz {dashed blue), and 150 GHz 
{dot-dashed red). 

Rather than parameterize a luminosity function, we di- 
rectly parameterize the mean comoving emissivity den- 
sity since that is all we need to make predictions of the 
clustered power. We write it as 

j{v,z) oc ax^(^;)cxp y 2a^) ■^''^^+^'>' ^^^^ 

In Figure |31 we show the contribution to the mean inten- 
sity as a function of rcdshift, dl^/dz — aj{v, z)dx/ dz, at 
150 and 1200 GHz (250 ^im) for both the LDP model and 
our fiducial singlc-SED model. We choose Td = 34 K, 
Zc ~ 2, az = 2, l3 = 2, and amid-iR = 2 as our fidu- 
cial single-SED model parameters for reasons given in 
Section 16.11 From Figure [31 we see that the redshift at 
which dlu/dz peaks decreases with increasing observing 
frequency, so the power seen by BLAST is predominantly 
coming from lower rcdshift than that seen by SPT. 

In Figure |4l we show how power spectrum shapes at 
150 GHz change with adjustments to Td and z^- We also 
show how they change with a switch from the single-SED 
model to the LDP model. All four curves are normalized 
to the same value si I = 3000. The relative percent dif- 
ference in power between singlc-SED model curves for 



various parameter values is small compared to the ac- 
curacy with which we can measure the power. Thus, 
we fix the template shapes at each frequency, calculated 
with our fiducial single-SED model, and allow only the 
amplitude to vary in this analysis. 

2.3. Clustered Radio Galaxies 

Models of galaxy evolution for both radio galaxies and 
DSFGs generically predict a clustering amplitude on the 
order of a few percent of the intensity of the mean back- 
ground on arcminute scales. Therefore, the amplitude of 
the mean background is the quantity of importance. The 
mean temperature due to undetected radio galaxies is 
difficult to estimate, and will be model dependent. How- 
ever, there are strong constraints from lower-frequency 
experiments. Detailed models by iGervasi et al.l (|2008f ) 
predict that the mean brightness temperature due to 
faint radio sources is 230mK(i^/GHz)~^'^ at low frequcn- 
cies, while recent re sults from the ARCADE2 experiment 
(jFixsen et al.ll2009l) find the mean extragalactic temper- 
ature above the CMB at frequencies from 3 to 10 GHz 
to be (1260 ± 90) mK (j//GHz)-2-6. The discrepancy be- 
tween source mode ls and the observed background is cur- 
rently unresolved (jSeiffert et al.ll2009D . but as an upper 
limit on the possible contamination from radio galaxies 
we can extrapolate the observed ARCADE2 temperature 
from their observed range of 3 — 10 GHz to 150 GHz, find- 
ing Tradio ^ 2.7 in antenna temperature, or about 
5 fiK in CMB temperature units. Assuming that 10% of 
this is clustered, the power contribution from clustered 
radio sources is about 0.25 fiK^. 

Based on small-scale CMB anisotropy measurements 
at 30 GHz from t he Sunyaev-Zel'dovich Array (SZA) 
(jSharp et al.ll2010l ) and the Cosmic Background Imager 
(CBI) (jSievers et al.|[2009l ). the above value of 0.25 fiK^ is 
almost certainly an overestimate of the contribution from 
clustered radio sources. In radio source models, the dom- 
inant contribution to the mean extragalactic background 
comes from faint sources, although there are substantial 
contributions from a broad range in flux. Extrapolating 
the ARCADE2 results to 30 GHz predicts the power from 
clustered radio sources to be of order 300 yuK^. This pre- 
diction is well above the SZA 95% confidence upper limit 
of 149 and would account for the entire measured 
CBI excess power. Assuming that the SZA upper limit is 
entirely due to clustered radio galaxies yields a clustered 
radio contribution of roughly 0.1 ^K^ at 150 GHz. Ex- 
trapolations from measured source counts predict a more 
reasonable contribution of ^ 6 /iK^ from clustered radio 
galaxies at 30 GHz (i.e., a factor of 50 lower than that in- 
ferred from the ARCADE2 extrapolation), which trans- 
lates into much less than 0.01 fj,K^ in power at 150 GHz. 
We can, therefore, safely neglect the clustering of radio 
sources. 

2.4. Unclustered Radio Galaxies 

Unlike DSFGs, the power contributed by radio sources 
is expected to be dominated by brighter sources; the 
number counts as a function of flux are such that the 
quan tity dl^/dlnS = S'^ dN/dS increases toward higher 
flux (|de Zotti et al.|[2005f ). The detected radio sources in 
the SPT maps are consistent with the expected distri- 
bution as shown in V09 and here in Figure [T] Applying 
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Fig. 3. — Redshift-dependence of contributions to the mean background intensity, dl^/dz at 150 and 1200 GHz (A = 250 fim) 
for the LDP model (blue solid curve) and our fiducial single^SED model [red dashed curve). For these figures, the single-SED 
model is normalized, at both frequencies, to have the same lu as the LDP model. 




1000 2000 3000 4000 5000 6000 

I 

Fig. 4. — Dependence on model parameters of the shapes of 
the angular power spectra of clustered point sources for our 
single-SED model at 150 GHz for various parameter combi- 
nations. We show three curves all with = 2, P — 2, and 
amid-iR = 2 and with the following Td and Zc combinations: 
= 34K and Zc = 2 {thm, solid, black line), T^j = 34K and 
Zc = 1 {dotted, blue line), Td — 50K and Zc = 2 {dashed, red 
line). Also shown for comparison is the shape of the LDP 
model {thick, solid, green line). All power spectra are nor- 
malized to the Td = 34K and Zc = 2 model at £ = 3000. 

Equation [2] to the Ide Zotti et al.l ()2005[ ) source counts, 
we calculate the expected power from radio sources be- 
low the 150 GHz SPT 5a flux threshold of 6.4 mJy, to 
be Ci = 7.9 X 10~^/zK^. We also examine the power 
spectrum of the all-sky radio source maps produced by 
S09 and find a very similar levels of residual radio source 
power {Ci = 8.0x IQ-'^/zK^) after masking sources above 
6.4 mJy. This is about 11% of the power we expect from 
DSFGs and is a significant source of uncertainty in our 
measurement of the mean spectral index of the DSFGs 
which we consider in Section 16.11 



2.5. Correlations with Galaxy Clusters 

A possible concern for results derived from the 150 GHz 
autospectrum is that radio galaxies and DSFGs are cor- 
related with galaxy clusters. In this case, the SZ effect 
from the galaxy clusters can "fill in" flux from the sources 
and lead to a suppressed signal. There are several lines 
of reasoning that argue against si gnificant radio c orre- 
lations; for example, the works bv iLin et al.l (|2009f ) and 
S09 show the number of cluster-correlated radio sources 
is expected to be small. We also test this assumption us- 
ing the simulated tSZ and radio source maps produced 
by S09 for 148 GHz. We compute the power spectra at 
i = 3000 for the tSZ map, the tSZ-|-radio source map, 
and the radio source only map after masking sources 
above 6.4 mJy. We find an anti-correlation coefficient 
of 2.3% for the two components which agrees well with 
other lines of argument. 

We also do not expect significant correlations between 
galaxy clusters and DSFGs. Galaxy clusters are not ob- 
served to be strongly overdense in DSFGs relative to the 
field (a fact or of 20, while t he mass is overdense by a fac- 
tor of 200 (jBai et al.ll2007D ). so only a small fraction of 
the galaxies contributing to the millimeter background 
can live in galaxy clusters^^. 

3. MODELING THE OTHER SIGNALS 

The clustered power and the Poisson power arc not 
the only significant contributors at 150 and 220 GHz for 
the angular scales of interest. In order to interpret the 
data, we also need to model fluctuation power from pri- 
mary CMB anisotropy, the tSZ effect, the kSZ effect, and 
Galactic cirrus. 

3.1. CMB 

■^■^ Note that this latter argument is inconsistent with the model- 
ing of DSFGs in S09. In their halo model, massive galaxy clusters 
would be expected to have roughly 200 times the density of sourees 
in the field, since the number of galaxies roughly tracks the mass 
of the halo. 
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We use CAMB (jLewis et al.ll2000f ) to predict the pri- 
mary CMB anisotropy for the standard, six-parameter, 
spatiaUy flat, lensed ACDM cosmological model. The 
six parameters are the baryon density fib, the density of 
cold dark matter £7^, the optical depth to recombination 
r, the angular size of the sound horizon on the surface of 
last scattering Q, the amplitude of the primordial density 
fluctuations Zn[10^''As], and the scalar spectral index Ug. 
We have not explored alternative cosmological models, 
although such models could affect the point source re- 
sults presented in this work by modifying the predicted 
high-£ primary CMB anisotropy. 

Gravitational lensing of primary CMB anisotropics by 
large-scale structure tends to increase the power at small 
angular scales which could impact the measurement of 
the clustered and Poisson point source contributions. In- 
cluding lensing out to ^ = 10, 000 proved computation- 
ally challenging. We avoid this computational constraint 
by calculating the lensing contribution for the best-fit 
CMB cosmological model, and adding this estimated 
lensing contribution to the unlensed CMB power spectra 
calculated for each parameter set. Given the small range 
of flm allowed by the combination of WMAP5, ACBAR, 
and QUaD data, we determine that using a fixed lensing 
contribution will misestimate the actual lensing by no 
more than 30%. The lensing contribution to the high-€ 
CMB spectrum is ~1.5 /iK^ near £ = 3000, and falls at 
higher i. This 0.5 ^K^ lensing error (30% of 1.5 ^iK^) is 
smaller than the clustered power we detect at 150 x 220 
and 220 x 220 by factors of 40 and 140 respectively. 

We also use EmuCMB^^ to calculate the lensed pri- 
mary CMB anisotropy. EmuCMB is a CMB Boltzmann 
code emulator that can return a lensed CMB tempera- 
ture power spectrum out to / = 5000 thousands of times 
faster than CAMB. We use it for the subtraction of the 
CMB contribution to the bandpowers in order to study 
residuals of our model fits, as described below. 

3.2. Thermal SZ 

The tSZ template we use is based on the simulated tSZ 
maps released by S09 for a WMAP5 cosmology. There is 
significant theoretical uncertainty in the shape and am- 
plitude of the tSZ effect for a given cosmology. The am- 
plitude uncertainty is mitigated by making the normal- 
ization of the tSZ template a free parameter at 150 GHz. 
We assume that the tSZ contribution to the 150 x 220 
and 220 x 220 bandpowers is zero which should be an 
excellent approximation at current significance levels. 

3.3. Kinetic SZ 

Galaxy clusters and other large-scale structures carry 
line-of-sight velocities that Doppler shift scattered CMB 
photons, resulting in a CMB te mperature anisotropy 
known as the kinetic SZ effect (jSunvaev fc ZeldovichI 
Il980f ). The signal is proportional to the product of 
the free electron density and velocity and has a spec- 
trum identical to a variation of the primary CMB black- 
body temperature. In contrast to the tSZ effect, elec- 
trons with temperatures as low as lO** K contribute to 
the kSZ effect. Therefore higher redshifts, before mas- 
sive objects finish collapsing, have a larger contribution 

•^^ Software and results of validation testing are available at 
[http:/ /www.emucmb.info^ 



to the kSZ effect. Clusters of galaxies at low redshift 
dominate the kSZ power on small angular scales, while 
high-redshift reionizing regions that are several tens of 
Mpc across have their largest relative contribution on 
angular scales around i ^ 2000, close to where sec- 
ondary anisotropics and point sources begin to dominate 
the anisotrop y power. This so-called patchy reionizat ion 
contribution |Gruzinov fc Hulll998t iKnox et al.lll998D is 
highly dependent on details of reionization and hence is 
highly uncertain. 

The kSZ contribution is an important source of un- 
certainty in our determination of power from clustered 
DSFGs. We first analyze the SPT bandpowers with a 
baseline assumption for the kSZ power, and then discuss 
the impact of uncertainty in the assumed kSZ signal on 
our conclusions. The baseline model is taken from the 
treatment of S09, and includes post-reionization contri- 
butions to kSZ calculated for WMAP5 best-fit cosmo- 
logical parameters. Patchy reionization may contribute 
significantly to the kSZ power spectrum and, therefore, 
we consider the implications of this potential additional 
kSZ power on our conclusions. 

3.4. Galactic Cirrus 

Another source of fluctuation power is emission from 
the cirrus-like dust clouds in our Galaxy. Although the 
area of sky used in this analysis has a relatively low level 
of emission from Galactic cirrus as m easured by the In- 
frared Astronomical Satellite (IRAS) (jNeugebauer et aP 
119841 ) , this emission is present at some level and its effect 
on the SPT spectra must be quantified. To this end, we 
cross-correlate the 150 GHz and 220 GHz SPT maps with 
predictions for cirrus emission at these frequencies from 
model eight of iFinkbeiner et al.l (|1999l hereafter FDS). 
We detect correlated power at high significance at both 
frequencies. The correlated power displays a Q oc £~^-^ 
shape after being corrected for the transfer functions of 
SPT and FDS. We use the normalization of the corre- 
lated power to predict the level of cirrus power in the 
SPT spectra. In practice this is implemented as a prior 
on the normalization of a oc £~^'^ term in the multi- 
component fit to the SPT power spectra. We find that 
the inclusion of this term has a negligible effect on all 
components except for the amplitude of the clustered 
point source power, which decreases by 8% (or O.Bcr) at 
150 X 220 and 7% (or 0.2ct) at 220 x 220. 

We note that the FDS predictions are based on 100 iim 
emission and extrapolated to millimeter wavelengths us- 
ing dust temperature estimates from the Diffuse Infrared 
Background Experiment (DIRBE) on the COBE satel- 
lite. DIRBE has degree-scale resolution, so if there is 
significant structure in the dust temperature at smaller 
scales, there will be some loss of correlation between the 
FDS template and the true cirrus emission in the SPT 
bands. This could bias our prediction for cirrus power 
at 150 and 220 GHz low. In order to estimate the mag- 
nitude of this effect, we have studied the angular power 
spectrum of the dust temperature in the FDS model and 
extrapolated it to smaller angular scales. We find that 
temperature structure could increase the predicted cir- 
rus power by up to 10%; such a small change in our prior 
on the cirrus power would have a negligible effect on our 
results, so we choose to ignore any decorrelation effects 
from temperature structure. 
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4. DATA COLLECTION AND DETERMINATION OF 
POWER SPECTRA 

We use data from one 100 deg^ field centered at right 
ascension 5^30"\ declination —55° (J2000) observed by 
SPT in the first half of the 2008 austral winter. We 
use the central 78 deg^ of the map which has uniform 
noise at a depth of 18 — arcmin at 150 GHz and 
40^K — arcmin at 220 GHz. Point sources detected at 
> 5cr (6.4 mJy) at 150 GHz arc masked from both maps 
in the analysis. 

From this two-frequency map, LIO estimate bandpow- 
ers for the 150 GHz auto spectrum, the 220 GHz auto 
spectrum, and the 150 x 220 GHz cross-spectrum. The 
bandpowers are estimated with a pseudo-C^ method. A 
cross-spectrum based analysis is used to eliminate noise 
bias ( fPolenta ct al .; 2005; Tristram et al. 2005j). Correc- 
tions for the time-stream filtering, beams and finite sky 
coverage are app lied using the fram ework of the MAS- 
TER formalism (jHivon et al.ll200^ . Bandpower uncer- 
tainties arc determined using signal-only simulations to 
estimate sample variance, and using the scatter between 
disjoint sets of cross-spectra to estimate the instrumen- 
tal noise variance. A full description of the bandpower 
estimation may be found in LIO and the bandpowers and 
window functions may be downloaded from the SPT Web 
site.24 



5. ANALYSIS OF THE THREE POWER SPECTRA 

We show the bandpowers for the three power spec- 
tra and maximum-likelihood models for the 150 x 220, 
and 220 x 220 power spectra in Figure [5] For 150 x 150 
the tSZ, kSZ, and clustered point source power spectra 
are "best-guess" models as explained below. The model 
includes the primary CMB, tSZ, kSZ, cirrus, and both 
Poisson and clustered DSFG components. In order to 
better visualize the point source contribution and need 
for a clustered component, we also plot the bandpow- 
ers after subtracting the estimated primary CMB, SZ, 
and cirrus contributions in Figure |B1 The bandpowers in 
Figs. [5] and [B] have been corrected for an ^-independent 
calibration factor and an ^-dependent beam factor (see 
LIO) for comparison with the models. Calibration and 
beam uncertainties have not been included in the error 
bars. 

Strictly for visual inspection (to guide the "x by eye" ) , 
the residual bandpower error bars in Figure [S] have been 
expanded to contain the uncertainty in the subtracted 
CMB contributions. This uncertainty is estimated by 
examining the scatter in a random sampling of CMB 
models from the parameter chains. Fixed tSZ (only at 
150 X 150), kSZ, and cirrus templates are also subtracted. 
This method is approximate in that we have ignored pos- 
sible correlations between these various components in 
doing the subtraction. 

From the center and right panels of Figure [Sj it is clear 
that we need the clustered component (or something like 
it) to avoid having an excess of power at low £ for the 
150 X 220 and 220 x 220 bandpowers. For a single spec- 
trum, the \ow-i excess could be explained by a sufficiently 
large kSZ effect. However, the relative powers seen in the 
150 X 150, 220 X 220, and 150x220 bandpowers are wefi 

■^'^ http:/ /pole. uchicago.edu/public/data/lueker09/ 



fit by a source with a dusty rather than a kSZ frequency 
dependence. 

Ideally, we would perform a simultaneous analysis of 
all three power spectra, which would help in discrimi- 
nation between the numerous components with similar 
spatial properties. It is clear from Figure [S] that such an 
analysis would give strong quantitative constraints. A 
multifrequency analysis with its additional challenges is 
planned for future work. In this work, we explore con- 
straints that we can set with the individual band powers. 
First, we consider constraints from the 150 x 150 band- 
powers alone. 

The expected contribution of tSZ and kSZ to our 
150 X 150 bandpowers complicates our analysis and ulti- 
mately prevents a definitive detection of clustered point 
source power at 150 x 150. The similarity of the as- 
sumed ^-space shapes of these three components leads to 
a near degeneracy in their amplitudes when constrained 
by the 150 x 150 bandpower data alone. Consequently, we 
perform an analysis with just the clustered point source 
power spectrum template, and treat it as a proxy for the 
power from all three contributions combined. We can 
interpret the constraints on the sum of these powers as 
an upper limit for any one of the signals. We find for the 
clustered power a.t £ = 3000 



1) 



27r 



-C 



PS,C 
150x150,; 



< 14 fiK^ 
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at 95% confidence. There are slight differences in the 
template shapes that lead to source-dependent upper 
limits. Using the S09 kSZ power spectrum as a template, 
we find the power at £ = 3000 to be < 13^X2.25 The 
analysis of LIO uses a DSFG-subtracted map and the as- 
sumption that the kSZ signal is given by S09 to produce 
the tighter upper limit for the kSZ power a.t £ = 3000 of 
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at 95% confidence. These estimates and upper limits 
are tabulated in Table [1] For these upper limits, we 
conservatively ignore the contribution from cirrus. 

In order to produce the upper limit on the clustered 
point source power at 150 x 150, we effectively assumed 
that the point source fiuctuations are not highly anti- 
correlated with total SZ. In the case of a large anti- 
correlation, it would be possible for the addition of 
sources to decrease the total fluctuation power. We 
might expect some degree of anti-correlation since more 
massive clusters will have increased (negative) tSZ flux 
and more member galaxies, some fraction of which will 
have dust and radio emission. As described in Sec- 
tion l2.5l the simulations presented by S09 have been used 
to constrain the correlation of tSZ with residual radio 
emission to a negligible level. LIO also argue that the 
anti-correlation of tSZ with dust emission is sufficiently 
small to have negligible impact on our upper limits. 

In addition to constraining the total power, we also 
make a "best-guess" at the amplitude of each of these 
three signals in the 150 x 150 data. We assume the kSZ 
model from S09. For the tSZ effect, we take the best-fit 

Our upper limits for both the kSZ shape and clustered point 
source power shapes are the same if specified at £ = 3800. 




Fig. 5. — The SPT bandpowers {black circles) are plotted at each frequency along with the following model components: 
total power {thm, solid, black line), CMB {thick, solid, purple line), tSZ eSect{dotted, blue line), clustered point source power 
{dot-dashed, red line), Poisson point source power {thin, dashed, cyan line), kSZ effect {thick, dashed, green line), and Galactic 
cirrus {thick, tripple- dot- dashed, orange line). In all cases, the kSZ effect is fixed to the S09 model. At 150 x 220 and 220 x 220, 
the plotted models are the maximum-likelihood solutions. We also plot the maximum-likelihood CMB and Poisson point source 
components for the 150 x 150 bandpowers; however, the other components plotted at 150 x 150 are extrapolated from other 
data. The tSZ component is set to the LIO maximum-likelihood amplitude for this kSZ model. The clustered point source 
component at 150 x 150 GHz is fixed to the value inferred by extrapolating from the best-fit amplitudes of the other two power 
spectra assuming perfect correlation. The cirrus component at 150 x 150 GHz is fixed to the value predicted by cross-correlating 
with FDS model eight. 
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Fig. 6. — The residual SPT bandpowers {black circles) after subtraction of the tSZ, kSZ, CMB, and cirrus model components 
shown in Figure [5] Note that the power is plotted in Ci oc Dt/{1{1 -\- 1)) as the Poisson point source power is flat in Ce. The 
residual model components are as follows (see the description in the Figure [5] caption) : total residual power {thin, solid, black 
line), clustered point source power {thick, dot-dashed, red line), and Poisson point source power {thin, dashed, cyan line). The 
error bars have been expanded to include uncertainties in the CMB contribution, as described in the text. 



tSZ amplitude for this kSZ model from LIO. The best- 
guess clustered point source power is an extrapolation 
from the best-fit clustered point source power values at 
220 x 220 and 150 x 220 assuming the S09 kSZ model 
and a correlation coefficient of unity between the spa- 
tial fluctuations due to clustering at 150 and 220 GHz. 
Wc assume the predicted level of cirrus discussed in Sec- 
tion l3.4l These estimates produce the tSZ, kSZ, clustered 
point source, and cirrus curves plotted in the left panel of 
Figure m Wc subtract these tSZ, kSZ, and cirrus models 
as well as the best-fit CMB model from the 150 x 150 



bandpowers to produce the left panel of Figure [6] 

6. COMPARISON WITH MODELS 

Here we make comparisons between modeled and ob- 
served frequency dependences for both Poisson and clus- 
tered power. We first examine the frequency dependence 
near the SPT observing frequencies. Our main conclu- 
sion is that the Poisson power favors models with steeper 
spectral indices. The steep frequency dependence of our 
detected Poisson power favors (3 = 2 over /3 = 1.5 and 
requires an adjustment to the spectral templates used in 
the LDP model. 
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TABLE 1 

Summary of SZ Results 



Signal 


Measured Vsooo (mK^) 


Upper limit on kSZ 


< 13 (95%) 


DSFG-subtracted 




tSZ + 0.46 kSZ 


4.2 (± 1.5) 


tsz 


3.3 (± 1.5) 


Residual clustering 


< 0.3 (95%) 


Fiducial kSZ 


2.0 



Note. — Upper limits and best-fit values evaluated at £ = 3000 for tSZ, kSZ, and clustered power at 150 GHz. The model 
templates for all three components are essentially completely degenerate at 150 GHz. The first row has the flat power determined 
from the 150 GHz bandpowers which consists of the sum of all three components. We report the upper limit on kSZ power at 
95% confidence from the 150 GHz maps in the second row, derived under the assumption that all fiat power is kSZ. The next 
three rows are results from the DSFG-subtracted maps in LIO which are listed here for comparison. The subtraction technique 
reduces the kSZ relative to the tSZ by a factor of ~ 2 and nearly eliminates clustered DSFG power as discussed in Section [T] 
The implied tSZ level under the assumption of the fiducial kSZ model is found on the fourth row. The upper limit at 95% 
confidence on residual clustered power in the DSFG-subtracted maps is on the fifth row. The fiducial kSZ model power is shown 
in the last row for reference. We conservatively ignore the expected cirrus contribution when estimating the kSZ upper limit. 
The quoted central values are the medians of each distribution. The quoted error bars are one half of the 68.3% confidence 
interval. 



TABLE 2 
Summary of Point Source Results 



Signal 


150 X 150 ifjK^) 


150 X 220 (/xK^) 


220 X 220 (AtK^) 


Q!150-220 


Poisson 
(Radio subtracted) 


r)Ps,p 

(jPS,P 


7.36 (± 0.59) xlO"'' 

10.5 (± 0.8) 
6.57 (± 0.70) xlO-" 


2.21 (± 0.16) xlO"^ 

31.7 (± 2.3) 
2.15 (± 0.17) xlO-^ 


7.88 (± 0.91) xlO"^ 

113 (± 13) 
7.84 (± 0.91) xlO"^ 


3.68 (± 0.20) 
3.86 (± 0.23 ) 


Clustered 


■'-'^3000 


< 14 (95%) 


16.4 (± 4.8) 


54 (± 18) 


3.8 (±1.3) 



Note. — Spectral indices and power for the Poisson and clustered point sources in the 150 x 150, 150 x 220, and 220 x 220 
spectra after masking the brightest sources. Recall that Vi = £{£ -\- l)Cf/(27r). The quoted central values are the medians of 
each distribution, and the quoted error bars are one half of the 68.3% confidence interval. The central values for the Poisson 
and clustered power are from chains with varying Poisson, CMB, clustered, and cirrus contributions and with the fixed fiducial 
kSZ (D3000 = 2.0 /iK^) contribution and the corresponding best-fit tSZ contribution from LIO (©3000 = 3.3 /iK^ at 150 x 150). 
For the Poisson term, af5o_220 is determined from all three spectra with an assumed upper limit on the variation of the spectral 
index of a a < 0.7 as discussed in the text. We also quote the effective power determination and spectral index fo r the Poisson 
component after subtracting an estimate of the residual radio source contribution using the Ide Zotti et all (pOOl ') model. The 
uncertainties are slightly larger in this case due to the assumed 50% uncertainty on the subtracted radio source power. The 
clustered power upper limit at 150 x 150 assumes no kSZ, tSZ, or cirrus contributions and is quoted at 95% confidence. The 
clustered component spectral index, ai5o-220i is determined from the clustered power determinations for the 150 x 220 and 
220 X 220 spectra with the assumption that the fluctuations are completely correlated between bands. 



Extending to higher frequencies, we find our single- 
SED model is able to reproduce the frequency depen- 
dence of the clustered power from our 150 and 220 GHz 
measurements to measurements at 1900 GHz. In fitting 
across this wide frequency range, we become highly sen- 
sitive to a degenerate combination of dust temperature 
and redshift distribution. We discuss the redshift distri- 
bution of the emissivity density used to match the data, 
and also compare to the shape of the mean spectrum 
inferred from COBE/FIRAS data. 

6.1. Spectral Indices at SPT Frequencies 

In Table [21 we show the amplitude of the Poisson point 
source power inferred from each of the three SPT power 
spectra. Assuming that each source has a spectral index 
between 150 and 220 GHz drawn from an independent 
normal distribution with mean Q!f5o_220 ^^'^ variance, 



tr^ , we can reconstruct the values of these two parameters 
from the Poisson amplitudes. The constraint on ct^ 
from the data is quite weak, and we apply a uniform prior 
of 0.2 < aa < 0.7 to determine af5o_220 = 3-68 ± 0.20. 
This prior on a a is motivated later in this section. 

This spectral index result does not include a correc- 
tion for the s ubdominant radio sou rce contribution. Ac- 
cording to the Ide Zotti et al.l (|2005l ) model, radio sources 
below our 5cr masking threshold will account for 11% of 
the reported Poisson power at 150 GHz and 0.5% of the 
Poisson power at 220 GHz. Subtracting the predicted 
radio power spectra with an assumed 50% uncertainty 

Note that 150 and 220 GHz arc only nominally the frequencies 
of these two (broad) channels. As described in the Appendix, the 
effective band centers, for the spectral shape expected for dusty 
galaxies, are 154.2 and 221.3 GHz; these are the frequencies we use 
for relating flux ratios to spectral indices. 
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steepens the measured spectral index and increases the 
uncertainty to give Q:f5o_220 ~ ^-^^ ^ 0.23. 

We also consider the impact of the uncertain shape of 
the "flat" power components on our inference of C^^'^ . 
We test the impact by fitting a single flat power compo- 
nent with a range of shapes together with a Poisson term. 
For the flat power component, we consider the range of 
tSZ and kSZ models in LIO in addition to the single-SED 
model in this work. The resultant, most likely C^^-^ 
values vary by 3 x 10^^/^K^ which leads to a spread in 
most likely '5i5o_220 "^f 0.07. This uncertainty has been 
added in quadrature to the (5f5o_220 values above and in 
Tabled 

Most models for millimeter-wavelength dust emission 
predict shallower spectral indices. For example, the Pois- 
son spectral index is 2.3 in the LDP model. The LDP 
model is in good agreement with published source counts 
at higher frequencies, and it is possible that the disagree- 
ment with the SPT data results entirely from the extrap- 
olation to our frequencies. All six of the template SEDs 
in the LDP model have ai5o-220 — 2.7 for z = 2. Ex- 
panding the set of spectral templates to include some 
which are steeper at long wavelengths might resolve the 
disagreement. 

We can also use our measurement of af5o_220 
constrain the dust emissivity index in the singlc-SED 
model. In the RayleighJeans (RJ) infinite-temperature 
limit a — 13 + 2, but finite temperature corrections to 
this relati on are significant. Ta king Td = 34 K (consis- 
tent with iChaoman et al.l [20051 who find = 36 ± 7 K 
and iDunne et al.l 120001 who find Td = 35.6 ± 4.9 K), 
we find ai5o_220 = /? + 1.7 for sources at 2; = 1 and 
ai50-220 = /3 + 1-5 for sources at z = 2. 

Of course, the emission is coming from a range of red- 
shifts, not just 2; = 1 or z = 2. If we take the red- 
shift distribution of the Poisson contribution to be the 
same as that of the mean background light in our fidu- 
cial single-SED model, we find af5o_220 = 1.38 -I- /3. This 
implies that /3 = 2.48 ± 0.23. Changing Td from 34 K 
to 40 K or changing Zc from Zc = 2 to Zc = 1, we find 
13 = 2.38 ± 0.23. However, BLAST and Spitzer data 
place some strong constraints on these variations as we 
describe below. 

Theoretical models calibrated with laboratory data 
and astronomic al observations at u ltraviolet and visi- 
ble wavelengths (jDraine fc LedlTOS^ lead to spe ctra well 
appro ximated at long wavelengths wi t h ~ 2 (jGordonI 
Il995| ). As reviewed in lDunne fc EalesI ()2001h . astrophys- 
ical observations of dust in different environments are 
consistent with /3 values between 1 and 2. Our data 
clearly prefer the high end of this range. We note that 
iMenv et"al1 (|2007[ ) find that effects due to long-range 
disorder in the dust grains can lead to /3 > 2 at long 
wavelengths. However, we emphasize that our results 
are consistent with (3 = 2 which is within 2.1cr of our 
best-fit value. 

Although (3 ^ 2is not surprising from the standpoint of 
theoretical calculations of the spectral properties of dust, 
many observations lead to constrain ts that appear incon- 
sistent with /3 = 2. For example, iDunne et al.l (j2000l ) 
find (3 = 1.3±0.2 for a sample of 104 galaxies selected at 
5 THz (60 ^ m) , and with fluxes m easured at 350 GHz and 
3 THz. And iKovacs et all (|2006[ ) find fr om observations 



of 15 galaxies that Td = 34.6 ± 3K(1.5//3)° '^i which for 
/3 = 2 is somewhat inconsistent with our assumed tem- 
perature of 34 K. However, relaxing the assumption of a 
single tempe rature can reconcile th e data with (3 = 2. 
For example, IDunne fc EalesI ()2001[ ) find a good fit with 
/3 = 2 by introducing a 20 K component to preserve the 
observed shallow spectral index from 350 GHz to 3 THz. 

We also compare our results with expectations from 
SEDs for which the spectral properties of the dust are 
calculated from first principles, and which includ e syn- 
chrotron and free-free emission (jSilva et al.lll998l here- 
after S98). The S98 templates were constructed from a 
detailed model of star and dust formation and calibrated 
with available data. In the left panel of Figure[7l we show 
the SEDs at z = 1 and z = 5 with luminosities set so that 
their fiuxes are 20 mJy at the observed v = 220 GHz. The 
single-SED model is also shown for comparison. In the 
right panel we show, as a function of redshift, the distri- 
bution of spectral indices derived from a uniform sam- 
pling of the S98 templates and reduced to a mean and 
la interval. The softening of spectral indices at lower 
redshifts is due to synchrotron and free-free emission. 
Assuming the distribution of template SEDs for galaxies 
is uniform, and that the redshift distribution of galaxies 
follows that of the 150 GHz background light of our fidu- 
cial single-SED model, we find a distribution of spectral 
indices with mean ai5o-220 = 3.04 and CTq = 0.28. This 
mean is inconsistent with our observationally determined 
value of Q!f5o_220 — 3. 86 ± 0.23. This difference is evi- 
dence that galactic spectra are preferentially more like 
the steepest S98 templates than the shallower ones. 

In order to constrain the residual DSFG Poisson power 
in the DSFG-subtracted map, LIO make use of a prior 
on the scatter in spectral indices of 0.2 < ctq < 0.7. Our 
exercise with the S98 templates would argue for ctq, ~ 0.3, 
but our difference with the S98 models argues for some 
caution and, therefore, we slightly more than double the 
dispersion for an upper bound of 0.7. Imposing a lower 
limit on (Tq is also conservative when used to place an 
upper limit on residual power; lower values of (Tq will 
produce less residual power. 

The DSFG clustering is subdominant to either the 
CMB anisotropy or Poisson DSFG power at all angular 
scales, and is therefore measured less precisely than the 
Poisson power. In addition, we have the difficulty of sep- 
arating the kSZ contribution which may have a similar 
shape. Table [2] reports clustered component amplitudes 
with the assumption that kSZ is at our fiducial value. As 
shown in the Appendix, the clustering spectral index of 
0^150-220 ~ ^-^ ^ ^-"^ reported in Table [5] is derived from 
the 150 X 220 and 220 x 220 bandpowcrs, and hence is 
independent of tSZ. However, the uncertain kSZ ampli- 
tude does introduce considerable extra uncertainty in the 
clustered DSFG amplitude and spectral index Q!i5o_220- 
If kSZ saturated our upper bound of 13 /iK^, we would 
find 

aj^Q_22o — 2.4. Assuming kSZ — 0, we find 

"f50-22o = 3-6± 1-2. 

We can also compare these spectral index results with 
those inferred for the DSFGs in SPT data that were re- 
solved out of the diffuse background light with S/N > 
4.5 at both 150 and 220 GHz and that do not have IRAS 
counterparts. Six of these sources are reported in V09. 
The posterior probability distribution of the spectral in- 
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Fig. 7. — Left hand side: The range of galaxy SED templates from ISilva et al.l (|1998l ) a.t z — 1 {blue shaded region) and z — 5 
{red shaded region), normalized to be 20mJy at 220 GHz. The SED for the fiducial single-SED model is also shown at 2 = 1 
{light blue line) and 2 = 5 {orange line). The shaded vertical regions show the positions and widths of the 150 and 220 GHz 
SPT bands use d in this work. Rig ht hand side: The 67% interval for 150 to 220 GHz model spectral indices at redshifts between 
and 6 for the ISilva et~aLl (|1998l ) galaxy SED templates. The fit to the Arp 220 SED from S98 {red, dot-dashed line) and the 
fiducial single-SED model SED {blue, dashed line) are also shown for comparison. 



dex for a source chosen at random from these six is ap- 
proximately fit by a Gaussian with 0:150-220 — 3.3 ± 0.7. 
The uncertainty includes noise and therefore represents 
an upper limit on the underlying scatter in the spectral 
indices of these sources. This is in reasonable agreement 
with the <5f5o_220 ^ 3.86 ± 0.23 calculated in this work 
for the fainter background population. 

Finally, we compare with observations made at higher 
frequencies but of objects with redshifts that are lower 
than the redshifts from which the bulk of the GIB light 
arises. These are of interest because they sample galaxies 
in the same rest-frame frequen cy range as probe d by the 
SPT observations of the GIB. iDve et al.l ()2009D present 
SED fits assuming (5 = 2 for 114 BLAST sources, mostly 
at z < 1, with 1.9 and 4.3 THz S pitzer observatio n s sup - 
plcmcnting the BLAST bands. iDunne fc Eale^ ()2001l ) 
find that the IRAS galaxy fluxes measured at 3 and 
5 THz by IRAS and at 350 and 670 GHz by SGUBA re- 
quire a two-temperature dust model with /3 between 1.5 
and 2. These observations which probe the same rest- 
frame frequency as the SPT observations, are roughly 
consistent with the high value of /? we deduce from the 
SPT data. 

6.2. Frequency- dependence of CIB Fluctuation Power 

In Figure [H we show the frequency dependence of 
Poisson power (in Jy^ sr~^) as inferred from SPT 
data at 150 a nd 220 GHz BLA ST data at 600, 860, 
and 1200 GHz (IViero et all I2009D . and Spitzer data at 
1900 GHz (jLagache et al.ll2007[) . In Figure H we do the 
same for the clustered power. The SPT data points are 
determined under the assumption of the fiducial kSZ am- 
plitude. Wc find that we can choose parameters for our 
singlc-SED model that reproduce the observed frequency 
dependence of the clustered power. 

Since the shapes of graybody spectra from individual 
objects at rcdshift z with temperature T are identical for 
equivalent T/ {1 -\- z), we have a degeneracy between the 
temperature and redshift distribution. Thus, we fix the 
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Fig. 8. — Amplitude of the Poisson point source power, 
C^^^^, as a function of i/cff = \/v\ x V2 inferred from SPT, 
BLAST (|Viero et al. 2009,), and Spitzer (|Lagache et al.ll2007D 
data. Also shown are the predictions from integrating the 
LDP model {squares) source counts up to the following flux 
cuts: 6.4 mjy at 150 GHz, 17mJy at 220 GHz (extrapola- 
tion of the 150 GHz flux cut with 0150-220 = 2.7), no flux 
cut at 600 GHz, 400 mJy at 860 GHz, 500 mJy at 1200 GHz, 
and 200 mJy at 1900 GHz. To compare the frequency depen- 
dence of the Poisson and clustered power, we also show the 
single-SED model clustered component frequency dependence 
{dashed line, this is equivalent to the shape of the solid line 
in Figure |9}. 

temperature when exploring the implications of our data 
in the context of the single-SED model. If the average 
temperature is higher than we assume, then the redshift 
distribution will be shifted toward lower redshifts than 
what we infer. 

For the parameter values from Section 16.11 (/3 = 2 and 
Td = 34 K), we find a good match to the data by setting 
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Fig. 9. — Amplitude of clustered point source power, 
, at £ = 3000, as a function of i/cS = V^^i x 1^2 inferred 



from SPT. BLAST jViero et al.|[2009l '). and Spitzer 1.9THz 
data (|Lagache et~aLl (20071 ). At 150 GHz the 95% confidence 
upper limit on the amplitude is given {down arrow). Also 
shown are predictions of the LDP model with bias set to 3 
(squares) and the same flux cuts as in Figure |8] Our single- 
SED model is shown {solid line) with amplitude and param- 
eters chosen to match SPT and BLAST power {Td = 34 K, 
Zc = 2, (Jz = 2, /3 = 2 and the following three values of 
Qmid-iR in order of decreasing p ower for comparison : 0.5, 
1, and 2). The mean inferred by IFixsen et al.l ()1998l ) from 
COBE-FYRAS data is also shown {dotted curve), which has 
a spectral dependence of Iv oc u'''^'^ B,y{18.5K) where Bi, is 
the Planck function. All inferences here (except for the upper 
limit) assume our fiducial value of kSZ is correct at 2.0 fiK^. 

Zc = 2 and 0-2=2. With these redshift-distribution 
parameters, the peak of the 150 GHz dl^/ dz is at z = 
3.1, and half the 150 GHz background light arises from 
z > 3.2. 

However, this inferred redshift distribution critically 
depends on the assumption of a single dust temperature 
of Td = 34 K. We know that the spectra of some galax- 
ies require significant additional cold components. Were 
we to include an additional (colder) dust component, we 
would find the redshift distribution shifted toward lower 
redshifts. Light from z — 2.4 and = 34K is spectrally 
the same as light from z = 1 and Td = 20 K. We do not 
include multiple dust temperatures in our modeling be- 
cause the small number of data points does not warrant 
such complexity. 

We can compare the dl^/dz of our model with that 
of the LDP model as a check on our consistency with 
known counts. We find that the mean redshifts of the 
GIB light in the LDP model and our single-SED model 
are quite similar: at 150, 600, and 1200 GHz the LDP 
model mean redshifts are 3.2, 2.2, and 1.6, respectively, 
and the fiducial single-SED model mean redshifts are 
3.3, 2.7, and 2.1. We can also compare our model 
with lMarsden et aLl (j2009D who find the fraction of back- 
ground light at z > 1.2 at 600, 860, and 1200 GHz to be 
61%, 51%, and 40%, respectively. This result relies on 
a stacking analysis using the FIDEL survey galaxies as 
a very deep catalog of sources of GIB fiux. Uncertain- 



ties in their result arise from possible incompleteness and 
from errors in the photometrically-determined division 
of FIDEL galaxies into z < 1.2 and z > 1.2 subcatalogs. 
In our fiducial single-SED model, there is significantly 
more light at z > 1.2 at these frequencies: 95%, 91%, 
and 85%, respectively. These differences may be due to 
the ab ove-mentioned uncertainties in the iMarsden et al.l 
(|2009[) result or they may be accommodated by a redshift 
distribution with more freedom than our two parameters 
allow, multiple SEDs, or multiple temperatures as dis- 
cussed above. 

Our single SED has the exponential decline on the 
Wien side of the graybody replaced by a power-law de- 
cline i/-"mid-iR as described and motivated in Section [2^ 
Adjusting amid-iR only affects the curve near the Spitzer 
data point in Figure [9l values Ofmid-iR = 0.5, 1, and 2 
(in order of decreasing power) are shown. 

As noted above, our single-SED model is not detailed 
enough to make predictions for Poisson power. However, 
we find that the observed Poisson power has a frequency 
dependence that is very similar to that of the clustered 
power. This similarity indicates that the spectral prop- 
erties do not strongly depend on galaxy brightness, red- 
shift, and bias. 

We can use our fiducial model above to relate Q;i5o_220 
to /3. Allowing higher Td or a distribution shifted toward 
lower redshift decreases /3 for fixed Q:i5o-220- However, 
the BLAST and Spitzer clustered power measurements 
limit our ability to relax the constraints on /3 in this 
way. By adjusting the temperature and the redshift dis- 
tribution, we can at best reduce (3 by 0.03 and maintain 
the same quality of fit. 

We find that a linear-theory bias of ~ 4 is required for 
the LDP model to ma t ch th e BLAST power, in agree- 
ment with IViero et al.l (|2009( ) . Such a bias would over- 
predict SPT power. The discrepancy could potentially 
be re solved by a frequency-dependent bias ()Knox et al.l 
I2OOH) . But the discrepancy between modeled and ob- 
served frequency dependence of the clustered power is 
similar to what we see with the Poisson power, and the 
Poisson power is unaffected by bias. Model adjustments 
to address the (more severe) Poisson power frequency- 
dependence discrepancy (such as the change to SEDs 
we suggest in this section) would, at the very least, 
go a long way toward eliminating the clustered power 
frequency-dependence discrepancy as well. Recall that 
the linear-theory bias is not a halo-model bias. Using 
a halo -model extension of the LDP model, IViero et al.l 
(l2009h find b ~ 2.4. 

IFixsen et all ()1998[ ) used COB£^FIRAS data to es- 
timate the mean of the baxikground light from 150 to 
2400 GHz. Their best fit is ^ = (1^/1^0)'^ 5^(18.5 ± 1.2 K) 
with /3 = 0.64 ±0.12. This form is to be interpreted phe- 
nomcnologically; i.e., it is not meant to represent emis- 
sion from an object with a single temperature. It can 
be understood as resulting from the sum of a number 
of graybody spectra with differing T/(l -I- z) and higher 
(3. The spread in T/{1 + z) values softens the rise in 
brightness with increasing frequency, since as frequency 
increases more objects are being sampled on the Wien 
side. The steeper rise we infer from SPT data does not 
contradict the FIRAS data which have no constraining 
power at 220 GHz and lower frequencies. 
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6.3. kSZ Models 

The contribution to kSZ from patchy reionization is 
dependent on the details of the transition from neutral 
to ionized inter-galactic medium and is thus highly un- 
certain. Our upper limit on kSZ is stringent enough to 
rule out the very highest predictions in t he literature. 



which range as high as 15 ^K^ at £ = 3000 ( Santos et al.l 



I2003f) ^'^. In more recent models fe.g.. iZahn et al.ll2005[ ). 
based on excursion set analytic models and full radia- 
tive transfer simulations, the patchy reionization signal 
is largest on scales of £ ~ 2000 where there is still con- 
siderable primary anisotropy power, and it is only tens 
of percents of the total kSZ signal on scales beyond the 
primaries. The modeling of the post-reionization kSZ 
has also been inconclusive, with amplitudes varying by 
a factor of 2. Given that we only marginally detect the 
thermal SZ effect which is generally expected to be larger 
than the kSZ signal, the current data quality is insuffi- 
cient to place a strong constraint on the amplitude and 
shape of the kSZ power spectrum. The increased sky 
coverage and addition of a 95 GHz channel in the 2009 
SPT data will make it possible to isolate and exploit this 
new cosmological probe. 

7. RESIDUAL CLUSTERED POWER IN THE 
"DSFG-SUBTRACTED MAP" 

Our 150 GHz maps have the following signal compo- 
nents 
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where PS, C stands for the clustered point source compo- 
nent and PS, P for the Poisson point source component. 
The 220 GHz maps are the same with the exception of 
no tSZ contribution. The DSFG-subtracted map of LIO 
is related to the individual frequency maps by 
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with X = 0.325; chosen to minimize the Poisson point 
source power. This value of x corresponds to 0150-220 = 
3.6 and is consistent with the spectral indices found in 
this work. A 1/(1 — cc) normalization factor can be in- 
cluded so that the subtraction has no effect on CMB 
temperature fluctuations, but we omit it here so that the 
subtraction has no effect on the tSZ signal. LIO analyze 
the power spectrum of this map, assuming it has negligi- 
ble contribution from clustered point sources, cirrus, and 
a constrained amount of residual Poisson power. Here, 
we justify the assumption of negligible clustered power. 

From the definition of the DSFG-subtracted map, the 
clustered power it contains is related to the auto and 
cross-clustered power via 
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(15) 

We do not measure the terms on the right-hand side well 
enough to use our measurements alone to provide a useful 
upper limit on the level of contamination. Fortunately, 
we have additional information to use. 

The power is for a model with r = 0.11; the predictions range 
even higher if we allow higher values of r. 



We assume that the spatial fluctuations at 150 GHz are 
100% correlated with the fluctuations at 220 GHz^^. For 
a power-law SED, the 150 and 220 GHz window func- 
tions would be identical, and we would have perfect cor- 
relation. Even for galaxies at z = 3, the SPT observ- 
ing bands are far on the RJ side of the spectrum and 
a power-law spectrum is a good approximation. There- 
fore, perfect co rrelation between bands is a very good 
approximation (|Knox et al.ll200l"l ). 

With this assumption, we have two estimates for the 
residual 150 GHz clustered power: 
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(17) 



where x^ is the true ratio between the 150 x 150 power 
and 220 x 220 clustered power, and x is the assumed 
value of 0.325. function of the spectral index 

0^150-220 of clustering component, which we assume 
(for the moment) to be equal to the spectral index of the 
Poisson term Q;f5o_220- approximate the distribution 
of spectral indices by a normal distribution with a = 3.86 
and CTq, = 0.23. 



We calculate the likelihood function for V 



PS.C 



from 



the 150 X 220 and 220 x 220 data. Assuming each estimate 
is independent, we combine the two likelihood functions 
to produce a combined estimate of the likelihood function 
according to 

P(I?ff''^) oc P(2?ff''^|150 x 220)P(X'ff''^|220 x 220) 
and integrate it to find the 95% confidence upper limit: 



PS.C 

,1 



< 0.3 /iK^ at ^ = 3000. This result has very httle 



V 

dependence on kSZ. This level of clustered point source 
power is small compared to the tSZ power. If we instead 
use the radio-unsubtracted spectral index of 3.68 ± 0.20 
to approximate the Xc distribution we get a lower 95% 
confidence upper limit of Pff'^ < 0.1 /iK^ at £ = 3000. 

We now consider how much cti^Q_22o might differ from 
Oif^Q_22o- We expect these two indices to be similar due 
to the fact that these are describing the spectra of light 
from the same objects. However, the light contributing 
to the determination of these indices comes from different 
weightings of these same objects with the Poisson power 
being more heavily weighted to brighter objects. One 
source of difference then is the dependence of the red- 
shift distribution on flux. To quantify how much different 
redshift distributions could affect the difference b e tween 
<^ 150-220 ^'^d af5Q_22o, we turn to the lRighi et al.l ()2008l ) 
model, which has very large amounts of the light com- 
ing from faint objects at high redshift. For this extreme 
model, the resulting clustered and Poisson indices shift 
with af5o_22o ~ "5f5Q_220 = 0.2. We thus conservatively 
test the effect of shifting the clustered power spectral in- 
dex with respect to the Poisson spectr al index by a s much 
as 0.4, double the difference in the iRighi et all ()2008l ) 
model. We repeat the upper limit calculation of the last 
paragraph with a mean spectral index of 3.46 instead of 



By which we mean (<5Ti5o<5T220> = \/ {^TI^q) {ST^2a) ■ 
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3.86 with the uncertainty unchanged. With this distribu- 
tion, we find a stronger upper hmit of 0.04 ^K^. Since we 
have very httle direct empirical knowledge of the sources 
in the faint magnitude ranges that make up the bulk of 
the CIB light at millimeter wavelengths, we adopt the 
original, more conservative assumption to quantify our 
uncertainty about residual CIB clustered power in the 
LIO DSFG-subtracted map. 

8. CONCLUSIONS 

We present detections of both Poisson and clustered 
power from DSFGs from an analysis of the SPT an- 
gular temperature anisotropy bandpowers presented in 
a companion paper (LIO). This is the first detection 
of the DSFG clustering component at millimeter wave- 
lengths. Wc find the Poisson power at ^ = 3000 to be 
10.5±0.8, 31.7±2.3, 113± 13 ^K^ for the 150 x 150 GHz, 
150 X 220 GHz, and 220 x 220 GHz power spectra, re- 
spectively. We find the clustered power amplitude at 
£ = 3000 to be approximately one half as large, with val- 
ues of 16.4±4.8 and 54±18 fiK'^ at 150 x 220 and 220 x 220 
respectively. At 150 GHz, where there is a degeneracy 
between the tSZ effect, kSZ effect, and clustered DSFG 
component, we set an upper limit on the clustered power 
of 14 /iK^ at 95% confidence. These probes of the statis- 
tical properties of the DSFG population complement the 
SPT source catalog presented by V09 which measures 
the bright tail of the DSFG and AGN fiux distributions. 

We also use the 150 GHz bandpowers to constrain the 
amplitude of the kSZ power spectrum to be < 13 /^K^ 
at £ = 3000 with 95% confidence. This upper limit be- 
gins to rule out the hi ghest of kSZ mode ls which reach 
15 nK"^ for T = 0.11 (jSantos et all [2001 . With more 
SPT data, including that from the 95 GHz band, we can 
expect tighter constraints on the kSZ effect and perhaps 
new information about the reionization history of the 
universe. 

We find the spectral dependence of the Poisson 
power due to DSFGs between 150 and 220 GHz to be 
Q;f5Q_22o = 3.86 ±0.23 which is on the steeper end of the 
expected range and higher than many current models. 
However, our results are consistent with theoretical ex- 
pectations for the spectral properties of interstellar dust 
and, under the assumption of no evolution of those spec- 
tral properties, with low-rcdshift observations of individ- 
ual galaxies at the relevant wavelengths. 

Combining the SPT results with those from BLAST 
and Spitzer, we show the frequency dependence of the 
Poisson and clustered power from 150 GHz to almost 
2000 GHz (160 /im). The clustered power over this large 
frequency range is consistent with a simple "single-SED" 
model of graybody spectra with a temperature of Td — 
34 K, (3 = 2, and a broad rcdshift distribution of the 
emissivity density, with half the light at 150 GHz com- 
ing from z > 3.2. We do not claim this model is a 
unique description of the current data, but over this 
broad range of frequencies the data tightly constrain the 
combination of dust temperature and redshift distribu- 
tion. The d u st tem perat ure of 34 K is c o nsiste nt with 
iDunne et all (|2000l ) and iChapman etHI ([20051 ). The 
redshift distribution is similar to that of the LDP model, 
but shifte d to higher red s hifts than the distribution in- 
ferred bv iMarsden et al.l ()2009D . We find it difficuh to 
accommodate values of /3 less than 2; reconciling this 



steep index with measurements of individual galaxies at 
shorter wavelengths may require dust with multiple tem- 
peratures. 

In a companion paper, SPT data are used by V09 to de- 
tect the brightest members of the source populations that 
are likely to contribute to these Poisson and clustered 
signals. Using the flux ratio between 150 and 220 GHz, 
these sources are separated into those with spectra con- 
sistent with active galactic nuclei (AGN) and those con- 
sistent with DSFGs. V09 present source counts for each 
population and argue that most of the detected DSFGs, 
having no counterparts in the IRAS catalog, must reside 
at high redshift or have unphysically cold dust. 

Furthermore, we argue that it is safe to neglect 
the clustering of radio sources given extrapolations to 
150 GHz of the clustered radio power from observa- 
tions of the mean temperature (suggesting less than 
^ 0.25 fJ,K^) and from source counts (suggesting less than 
0.01 /iK^). Several lines of reasoning also suggest that 
the correlations of radio galaxies and DSFGs with the 
SZ signal due to galaxy clusters will be negligible. 

LIO derive constraints on SZ power from a DSFG- 
subtracted map constructed from a linear combination of 
the 150 and 220 GHz maps. A residual clustered DSFG 
component would bias the SZ constraints presented by 
LIO. The residual DSFG power can be predicted from 
the average spectral index of the DSFG population, the 
measured clustered power amplitudes at 150 x 220 and 
220 X 220, and the exact frequency combination used to 
construct the DSFG-subtracted bandpowers. We assume 
that the frequency dependence of the clustered compo- 
nent is the same as the Poisson component. With this as- 
sumption, we place an upper limit at 95% confidence on 
the clustered component in the DSFG-subtracted band 
powers of 0.3 ^K^ at i = 3000. We also argue that the 
variation of spectral indices between sources is expected 
to be small; ctq G [0.2,0.7]. This constraint on spectral 
index variability is converted into a prior on the residual 
Poisson DSFG component in the DSFG-subtracted maps 
in LIO. 

A full multifrequcncy analysis is a natural progression 
from the single-frequency bandpower parameter fits per- 
formed in this work, and we are in the process of devel- 
oping a simultaneous multifrequcncy fitting pipeline. In 
addition, the SPT bandpowers used in this work repre- 
sent only a small fraction of the complete SPT data set. 
The complete SPT survey will cover over 2000 deg^ at 95, 
150, and 220 GHz. The three frequencies will allow us to 
separate the kSZ from DSFGs and tSZ and place new 
constraints on the reionization history of the universe. 
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APPENDIX 

RECONSTRUCTION OF SPECTRAL PROPERTIES 

The spectra used in this analysis are caUbrated in CMB temperature units. We can relate temperatures to flux 
according to 

5Iu = ^\tcms5T.- (A1) 

If we assume that the sources each have a spectral index a drawn from an independent normal distribution with 
mean a and variance cr^ then 

5Iu = 5h, ( — 



[l + 5a\n{v/v,)]. (A2) 



Taking v = vi^q and v^, = we determine the Poisson power spectra from two maps: 

{6ll^,) = {51^) 

V ^^220 / 

{5I150SI220) = (<5/')(i^l50/i'220)"- (A3) 
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We can directly solve this system to get the three unknowns in terms of the three known quantities: 



^150-220 - 



220/ 

ln((,5/i50(5/220)/(^/|20)) 

ln(i'i5o/'^22o) 



cr„ = [ln(j^i5o/;^22o)]" 



- 1 



(A4) 



In practice, we find a significant degeneracy between 01150-220 ^nd cr^ in the current data and choose to place a prior 
on CT^ for the q;i5o-220 constraints in this work. 

The spectral index variation only contributes to the Poisson component, not the clustered component^^. To param- 
eterize the frequency dependence of the clustered component, we introduce an effective spectral index ai5Q_220 ^^"^ ^ 
cross-correlation coefficient r, such that 

/ 7^1 rn \ ^"150-220 

('5/|2o) = («?5o) — (A5) 

V 2^220 / 

{SI150SI220) =rsJ{Sl!,o){Slho)- (A6) 

The value of r depends on the overlap of the window functions. For the LDP model, r ~ 0.998 at £ = 3000, with a 
very weak ^-dependence. 

In LIO; the DSFG-subtracted map is proportional to STi^q — x5T22o- If this linear combination is to minimize the 
contribution from a component with spectral index Q!i5o-220 then 

"^V<^^-o)/('5ri5o)-^^^(^^j =l-2l(^^j . (A7) 

The SPT bands only nominally measure power at 150 GHz and 220 GHz, and instead have a width of about 32% 
and 22%, respectively. For purposes of estimating spectral indices, we define effective band centers {ves) such that 
a power-law flux {S oc v") integrated over each band has a ratio S1/S2 ~ {i^cs .1 / ^cS ,2)°' ■ These band centers are 
weakly dependent on a, but differ significantly from the nominal values of 150 and 220 GHz. The effective band 
centers for a = 3.5 {a = 0) are 154.2 (152.0) GHz and 221.3 (219.5) GHz. In Section i] and in the Appendix, we 
use 1^150 = 154.2 GHz and 1^220 = 221.3 GHz to relate ratios of the 150 GHz band to 220 GHz band fluxes to spectral 
indices. Similarly, for the dB/dT factors in Equation IA7I we use the band-averaged values. 



This is true if we assume the spectral index scatter is uncorrelated with the large-scale number-density fluctuations. 



